A hydrodynamical trigger mechanism for pulsar glitches 
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The standard explanation for large pulsar glitches involves transfer of angular momentum from an 
internal superfluid component to the star's crust. This model requires an instability to trigger the 
sudden unpinning of the vortices by means of which the superfluid rotates. This Letter describes a 
new instability that may play this role. The instability, which is associated with the inertial r-modes 
of a superfluid neutron star, sets in once the rotational lag in the system reaches a critical level. 
We demonstrate that our model is in good agreement with observed glitch data, suggesting that the 
superfluid r-mode instability may indeed be the mechanism that triggers large pulsar glitches. 



Introduction. — Even though pulsars are generally very 
stable rotators, in some cases with an accuracy that rivals 
the best terrestrial atomic clocks, many systems exhibit 
a variety of timing "noise" . The most enigmatic features 
are associated with the so-called glitches, sudden spin- 
up events followed by a relaxation towards steady long- 
term spin-down. Several hundred glitches, with magni- 
tude in the range Ail c /il c w 10~ 9 — 10~ 6 where f2 c is the 
observed rotation frequency, have now been observed in 
over 100 pulsars The archetypal glitching pulsar is 
Vela, which exhibits regular large glitches. Glitches have 
also been reported in several magnetars as well as one, 
very slowly rotating, accreting neutron star Q. 

Despite the relative wealth of observational data our 
theoretical understanding of glitches has not advanced 
considerably in recent years. The standard "model" for 
large pulsar glitches envisages a sudden transfer of angu- 
lar momentum from a superfluid component to the rest of 
the star Q , which includes the crust (to which the pulsar 
mechanism is assumed rigidly attached) and the charged 
matter in the core. A superfluid rotates by forming a 
dense array of vortices, and the vortex configuration de- 
termines the global rotation. The key idea for explaining 
glitches is that, if the superfluid vortices are pinned to 
the other component, a rotational lag builds up as the 
crust spins down due to electromagnetic braking. Once 
the rotational lag reaches some critical level, the pinning 
breaks. This allows the vortices to move, which leads to 
a transfer of angular momentum between the two com- 
ponents and the observed spin-up of the crust. 

Most theoretical work has focused on either the 
strength of the vortex pinning 0, Q or the post-glitch 
evolution 7]. There have not been many suggestions for 
the mechanism that triggers the glitch in the first place. 
It is generally expected that this role will be played by 
some kind of instability, but there are few truly quantita- 
tive models. The results presented in this Letter change 
the situation dramatically. We present evidence for a new 
instability, acting on the inertial modes of a rotating su- 
perfluid star, that sets in beyond a critical rotational lag. 
The predictions of this model agree well with the obser- 



vational data making it plausible that this mechanism 
provides a missing piece in the pulsar glitch puzzle. 

Inertial mode analysis. — We want to improve our un- 
derstanding of the global hydrodynamics associated with 
a pulsar glitch. Even though this should be a key issue, it 
has not been discussed in detail previously. In principle, 
one would expect to be able to express the dynamics in 
terms of global oscillation modes of the system. In this 
Letter we present the first results for inertial modes of 
neutron star models with the two main features required 
for the standard glitch models, a superfluid neutron com- 
ponent that rotates at a rate different from that of the 
crust and pinned neutron vortices. 

We use the standard two-fluid model for superfluid 
neutron stars (see for example 0]), identifying the two 
components with the neutron superfluid and a conglom- 
erate of all charged particles (the "protons"). In the 
following, the index x = {n, p} identifies the distinct 
fluids. Our aim is to model small amplitude oscilla- 
tions with respect to a background configuration where 
both fluids rotate rigidly with (parallel) angular velocities 
'x = £lJ ^j x k and where the magnitudes are different, 
7^ ^p- The linear perturbations of this system (as- 
suming a time dependence ~ exp(io~t)) are, in the inertial 
frame, described by the two coupled Eulcr equations, 

(ia + viVj) 5v l n + <KVX + V*<ty» = SfLf (1) 



(ia + ujVj-) Sv p + SvlVjV z p + V*#p = 



(2) 



Here Sip x = Sjl x + SQ represents the sum of the perturbed 
specific chemical potential and gravitational potential. 
We have also introduced x p = p p j p n . This ratio, which is 
roughly equal to the proton fraction, is assumed constant 
throughout the star. For simplicity, we assume that the 
two fluids are incompressible, which means that Vidv^ = 
0. In general, the two fluids are coupled i) chemically, 
ii) gravitationally, iii) via the entrainment effect, and iv) 
by the vortex mediated mutual friction /' nf . For clarity, 
we will ignore the entrainment in the present analysis. 
A detailed discussion of how the entrainment affects our 
results will be provided elsewhere. 
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For the inertial modes, the main coupling mechanism 
is provided by the mutual friction force. The general 
expression for this force is 



B'e ijk cjiw* (3) 



where 



v p and uj 1 „ 



= e ijk WjVl. A "hat" de- 
notes a unit vector. When the two fluids are not coro- 
tating, the perturbed force 5f^ f is quite complex [loj ]. 
The form J3]) for results from balancing the Mag- 
nus force that acts on the neutron vortices and a re- 
sistive "drag" force which represents the interaction be- 
tween the vortices and the charged fluid. Representing 
the drag force by a dimensionless coefficient TZ, one finds 
that B' = TZB = TZ 2 /{1 + TZ 2 ). The range of values 
that TZ takes in a neutron star is not well known. The 
standard assumption in studies of neutron star oscilla- 
tions has been that the drag is weak, which means that 
8'c6< 1 . Then the second term in ([3]) has no effect on 
the dynamics. However, it may well be the opposite limit 
that applies. The vortices in a neutron star core may ex- 
perience a strong drag force if their interaction with the 
magnetic fluxtubes is efficient 11, 12]. The drag on the 



superfluid in the crust may also be strong due to vortex 
"pinning" by the lattice nuclei. Even though the current 
evidence [H, [y] favours weak crustal pinning, the existence 
of strong pinning regions has not been completely ruled 
out. In these cases one must consider the strong coupling 
limit TZ > 1 which leads to B « and B' ps 1. 

The hydrodynamical equations CE}-© allow for a rich 
set of oscillation modes. Our focus will be on a subset 
of the inertial modes, the purely axial r-modes. The r- 
modes have attracted considerable attention since they 
may suffer a gravitational-driven instability (l3| . They 
have been studied in superfluid neutron stars previously 

H 
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la], but not accounting for both mutual friction 
and a rotational lag. When expressed in terms of spher- 
ical harmonics, the r-mode velocity fields take the form 



= JmUfY? 



r 2 sin 2 i 



Ufd e Y{ 



(4) 



In the corresponding single fluid problem such purely ax- 
ial solutions exist, corresponding to a single I — m multi- 
pole. The solution (0| satisfies the continuity equations 
automatically, which means that we only have to con- 
sider the two Euler equations. Inserting (|4]) in (JTJ)- {2j) and 
eliminating Sip x one arrives at two equations for Uf. As 
in the single fluid inertial mode problem, these relations 
feature couplings between different / multipoles. In the 
present problem these couplings arc further complicated 
by the presence of Sf^ f . As a result, one would expect 
the problem to be difficult to solve. Remarkably, this is 
not the case. It turns out that, even when Q n ^ J7 p , there 
exists a simple I — to solution of the form = A x r m+1 . 
This solution is exact at leading order in rotation and for 
all values of TZ. The amplitudes A x follow from the solu- 



tion of a 2 x 2 algebraic system, the determinant of which 
provides the dispersion relation for the mode frequency. 

A new superfluid instability. — The strong coupling 
limit, B = and B' = 1, provides a good illustration 
of the main new result. In this case we obtain a simple 
mode solution. In terms of the dimensionless parameters 



if2 r 



A = 



we have the two frequency solutions 



1 



«1,2 



(to + l)x p 



(5) 



(G) 



where 



V = (l+a; p ) 2 +2A [1 + x p {3 - to(to + 1)}]+0(A 2 ) (7) 
The amplitudes are related according to 



A n 2 (1 + A) 



A v 



(to + 1)k 



(8) 



Let us focus on the short-lengthscale modes. Tak- 
ing to 3> 1 and recalling that x p and A are both small 
(generally A <C x p ), it is easy to show that one of the 
above r-mode solutions becomes unstable (Idi[k] < 0) for 
to > to c , where 



1 



y/2x~A 



320 



0.05 



1/2 



10 



-4\ V2 



(9) 



For m 3> to c the instability growth timescale Tg row 
l/(f2 p Im[K]) is well approximated by 



' grow 




— Wtt-t- ~ 0.25 



X n \V2 no 



0.05 



1/2 



0.1 s, 
(10) 

where P — 27r/ Sl p is the observed spin period (we asso- 
ciate the charged component with the crust, Sl p = fi c ). 
We see that this new instability can grow rapidly, on a 
timescale comparable to the rotation period of the star. 

Although we cannot yet claim to understand the de- 
tailed nature of this new r-mode instability, we have some 
useful clues. Firstly, we find from ((8]) that the unstable 
modes are such that |^4 n /A p | ~ x p . Thus, the fluid mo- 
tion is predominantly in the proton fluid. There should 
also be a close connection with the short wavelength in- 
stability that we recently discovered for precessing super- 
fluid stars [13] • In that case, the result followed from a lo- 
cal plane- wave analysis of the inertial modes. An attempt 
to link these two results would be useful. Finally, since 
the present system has two distinct rotation rates one 
might expect the instability to belong to the general two- 
stream class (l6| . Such instabilities are generic in multi- 
fluid systems. The intuitive condition for such an insta- 
bility dictates that the mode's pattern speed — Re(cr)/m 
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should lie between f2 n and Q p [l6[ ■ This criterion trans- 
lates into m > s/2m c , which is satisfied in the main part 
of the parameter space where the instability is operative. 

In order for the new instability to affect the dynamics 
of realistic neutron stars it must (at least) overcome vis- 
cous damping. For young and mature neutron stars, dis- 
sipation is dominated by shear viscosity due to electron- 
electron collisions [18j |. For a uniform density star with 
M = IAMq and R = 10 km (the canonical values used 
in the following) the corresponding viscosity coefficient is 






where T is the core temperature. We can use the stan- 
dard energy-integral approach to estimate the viscous 
damping timescale. In fact, since \A n \ <C |yl p | we can 
use existing results for r-modes of uniform density stars 
pol ]. remembering that shear viscosity only acts on the 
proton fluid. Thus, a simple calculation leads to 

3Mx p 6 x 10 4 /0.05\ 1/2 / T \ 2 
Tsv ~ 8irm 2 r) cc R ~ m 2 \~ J [Wk J S 

(12) 

Combining (Ti"Tj|) and (TT2")) we have a criterion for the un- 
stable modes to grow fast enough to overcome viscous 
damping. The condition r grow < r sv leads to 

/0.05\ 1/2 / A \ 1/4 / T \ /0.1s\ 1/2 

m < 500 (-j (i5=0 {wk){— ) 

(13) 

This gives the range of unstable m modes. Since the 
instability sets in when m > m c we conclude that the 
system becomes unstable once it reaches the critical lag 

/ p \ 2 / 3 / rp \ -4/3 

A ^ 6xl „->(_) (j^) ,14) 

Making contact with observations. — Within a two- 
component model, it is straightforward to estimate the 
critical lag required to explain the observations. Assum- 
ing that angular momentum is conserved in the process, 
one must have I c Afl c w — 7 s Af2 s , where I s and I c are 
the two moments of inertia, while AQ S and Af2 c rep- 
resent the changes in the corresponding spin frequen- 
cies. The glitch data suggests that about 2 % of the total 
spin-down is reversed in the glitches [2l[ , suggesting that 
I s /I c ss 0.02. In order to permit Vela-sized glitches with 
Af2 c /i7 c ~ 10~ 6 we then need (assuming fl c — f2 s after 
the event) 

A 9 «^«5xl0- 4 (15) 

The observational estimate of the lag A g at which large 
glitches occur is close to our estimate (fT4")l for the onset 



of the superfluid r-mode instability. This is unlikely to 
be a coincidence. Even though it is difficult to compare 
the parameters of our two-fluid neutron star model to 
the global quantities used in the phenomenological dis- 
cussion directly, it is clear that our new instability has 
the features expected of a glitch trigger mechanism. It 
operates in the strong drag limit, where vortices are ef- 
fectively pinned to the charged component. As long as 
the system is stable, a rotational lag should build up as 
the crust spins down. Once the system evolves beyond 
the critical level (Tj"4")) a range of unstable r-modes grow on 
a timescale of a few rotation periods. We cannot yet say 
what happens when these modes reach large (non-linear) 
amplitudes, but it seems inevitable that the fluid mo- 
tion associated with the instability will break the vortex 
pinning, allowing a glitch to proceed. 

Let us compare the "predictions" of our model to the 
data for pulsars exhibiting large glitches. To do this, we 
estimate the maximum glitch size allowed if A g = A c , 
assuming a completely relaxed system and I s /I c = 0.02. 
Since we do not have temperature data for most glitching 
pulsars we estimate T by combining the heat blanket 
model from with a simple modified URCA cooling 
law. Calibrating this model to the Vela pulsar, for which 
T w 6.9 x 10 7 K [H, we find 

/ t \~ 1/6 
T w 3.3 x 10 8 ( -j-^ J K (16) 

Here t c — P/2P is the characteristic pulsar age. The 
results are shown in Fig. [T] This Figure shows that our 
model does well in predicting the maximum glitches one 
should expect. The data is consistent with the idea that 
a system needs to evolve into the instability region be- 
fore a large glitch happens. It should be noted that, 
even though the instability first appears at m — m Cl the 
growth time is much longer than (|10[) until m > 1.2m c 
or so. It is also interesting to note that two of the sys- 
tems with actual temperature data [23|], Vela and PSR 
B1706-44, both sit on the m « 1.6m c curve. 

Discussion. - We have described a new instability 
that may operate in rotating superfluid neutron stars. 
We have demonstrated that this instability sets in at 
parameter values that compare well with those inferred 
from pulsar glitches. This suggests that the superfluid 
r-mode instability may be the mechanism that triggers 
large pulsar glitches. 

The model is consistent with a number of observed 
properties of glitching pulsars: 

i) Adolescent pulsars, like Crab and PSR J0537-69, 
should only exhibit small amplitude glitches. For fast 
spin and a relatively high temperature the instability sets 
in at smaller values of A. 

ii) More mature and slower spinning neutron stars 
have colder cores which means they can produce larger 
glitches, provided the required A can build up (cf. PSR 



4 



u 

< 



10" 



10" 



10" 

Crab 

7 



10 





1 ' 


2.5m . ; 




B2334+61____ 


7.6m 






m 

— — 


: UrO!. :.|..- A 




i% ; 


/Vela**' A 


A 


A 0.5% - 


• •Fi { 

7//v>'a .. 


A 




I'*' - 


A 
A* 

1 


a *; 

A - 







10 



(P/O.ls) (T/10 K) 



FIG. 1: The maximum glitches predicted from (|14[) are com- 
pared to observations. The m c curve represents the onset of 
the superfluid instability (the instability region is grey, and for 
2.5m c there is a range of unstable inertial modes) and a glitch 
involving the relaxation of 2% of the total moment of inertia 
(corresponding curves for 1% and 0.5% are also shown). The 
circles are systems with temperature data (taken from [23l]). 
The triangles represent systems for which T is estimated using 
(fT6l) . The glitch-data is taken from P,H^]. 



eters (P ~ 10 s, T ~ 10 9 K), we would not expect these 
objects to exhibit large glitches. Yet, they do [2]. The 
resolution may be that these glitches involve a larger frac- 
tion I s /I c - A very interesting question for the future con- 
cerns whether this is a natural consequence of stronger 
magnetic pinning in the neutron star core. 

These first results are promising, but we are obvi- 
ously far away from a complete understanding of this 
new mechanism. Future work needs to consider more 
detailed neutron star models. We need to understand 
the local mutual friction parameters and the nature of 
vortex pinning. The inertial mode problem for realistic 
neutron star models presents a real challenge. We should 
also make more detailed attempts at understanding the 
observations. A natural first step would be to improve 
on the temperature estimates used in Figure [TJ Finally, 
we need to understand the nonlinear development of the 
new instability. This problem can perhaps be studied 
with numerical simulations, building on the work dis- 
cussed in 26j. Ultimately one would hope to arrive at a 
truly quantitative model for pulsar glitches. 



J1806-21 with An c /n c ps 1.6 x 1CT 5 [24]). Since A c in- 
creases as the star ages one would expect neutron stars 
to cease to glitch eventually. 

iii) For any glitch mechanism that relies on a critical spin 
lag between a superfluid component and the rest of the 
star it is easy to estimate the time interval t g between 
successive glitches. Assuming that each glitch relaxes 



the system, we estimate t g 



2 ^ g f p*. t 2 ^ f ~fj f 



For Vela 



we then find t g > 750 d which compares well with the ob- 
served averaged time of about 1000 d. The most regular 
known glitcher, PSR J0537-69, has an average interglitch 
time of about 120 d In the absence of temperature 
data we use (jTHJ) for this object and find t g > 90 d. Again, 
the agreement with the observations is good, and consis- 
tent with the notion that the system evolves into the 
unstable regime before a glitch occurs, 
iv) There is no reason why the instability should not oper- 
ate in all spinning neutron stars in which a rotational lag 
builds up. In particular, one may expect accreting neu- 
tron stars to "glitch" occasionally. So far, there has only 
been one suggested event, in the slowly rotating transient 
KS 1947+300 S\. For this system our model suggests 
(combining P = 18.7 s with T ps 10 s K, a temperature 
that should be typical of an accreting star) a maximum 
glitch of Afl c /Q c ps 4 x 10~ 5 . This is very close to the 
suggested observed glitch with AQ C /Q C ps 3.7 x 10 -5 
Such events should, of course, be extremely rare. 

Since we have considered the non-magnetic inertial 
mode problem, our model does not apply (without mod- 
ification) to magnetars. It is nevertheless interesting to 
consider these systems. Given typical magnetar param- 
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